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ABSTRACT 


Atmospheric density tf LUCLUATI One induce phase 
perturbations that degrade the spatial coherence of a laser 
beam propagating through the atmosphere. These 
degradations spread the laser beam and alter the centroid 
and intensity profile stochastically. Turbulent conditions 
arising from various environmental Situations and 
meteorological phenomena are found at virtually all levels 
of the atmosphere. A substantial fraction of the optical 
turbulence along a vertical path arises from the surface 
heat flux within the first 100-200 meters above the ground. 

This thesis seeks to measure and analyze these 
turbulent layers. A high frequency acoustic echosounder 
was developed to analyze atmospheric turbulence within the 
first 100-200 meters above the ground with extremely high 
accuracy. The echosounder design incorporated a 25 
element, square, planar array housed within a lead and foam 
insulated, acoustic shroud. This acoustic echosounder was 
employed to collect real-time, low level atmospheric data 
at two sites. This data should prove useful in determining 
whether the performance of laser and electro-optical 
systems could be increased by raising sensors) and 
transmitters well above the local turbulent layer through 


the use of towers or other such structures. 
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I. INTRODUCTION 


A coherent laser beam propagating through the 
atmosphere 1S very susceptible to numerous turbulent- 
dependent processes. As electromagnetic waves transit the 
turbulent regions, atmospheric irregularities randomize the 
amplitude and phase of the wave. A high resolution 
acoustic profiler was designed and developed to quantify 
the altitude dependence of these atmospheric disturbances. 
Acoustic energy is more susceptible to these irregularities 
than electromagnetic waves. Historically acoustic sounders 
have proven to be the best means of providing high 
resolution atmospheric turbulence profiles within the first 
500 to 1000 meters above the ground. Acoustic echosounders 
are frequently used to detect and measure atmospheric 
density and velocity fluctuations resulting from wind 
shears, convection, and temperature inversions. 

Presently two optical parameters, the spatial coherence 
length (ro) and the isoplanatic angle (80) are measures of 
the perturbation of electromagnetic waves propagating 
through the atmosphere, and are measured very reliably by 
optical systems developed by Walters [Refs. Pana 2 |. 
Although these systems measure a line integral of the 
atmospheric optical turbulence with high accuracy, a major 


drawback of each system is that no provision is provided 


for determing the height of the atmospheric disruptions. 
If these disturbances are found to exist very near the 
surtace, it may be possible to negate their etfect by 
elevating sensors and transmitters or altering the surface 
heat flux to suppress the generation of turbulence. 

This thesis deals with the design, construction and 
implementation of a high frequency, acoustic echosounder 
which will accurately analyze the atmospheric density 


fluctuations within approximately the first 200 meters of 


the atmosphere. As this project is a product of the 
research and efforts of two students, the work was 
appropriately divided. My particular task was concentrated 


in the actual design and hardware development of the 
echosounder. My colleage, LT. Michael Wroblewski [{Ref. 3], 
dedicated his etforts toward software development and data 
display, analysis and storage. 

The echosounder developed and discussed herein utilizes 
a 25 speaker, square, planar antenna array enclosed within 
an acoustically absorbent shroud comprised of lead and foam 
insulation inside a sturdy plywood framework. This antenna 
array was employed as both the transmitter and receiver of 
the acoustic signal. A computer controlled function 
generator transmits an acoustic pulse through the antenna 
array. Additional electronic components incorporated 
within the overall echosounder system are used to amplify, 


digitize, display, analyze and ultimately store the return 


Signal. Over time the received signal data is used to 
produce a quantitative, three dimensional map of the 
atmospheric turbulence. 

Acoustic echosounders have been developed oy in use 
for many years and have proven to be a valuable probe for 
analyzing the structure and dynamics of the lower 
atmosphere [Refs. 4 through 7]. Devices similar to ours 
have been used to obtain profiles of the atmospheric 
density and temperature fluctuations [Refs. 8 and 9]. Our 
device excels by utilizing a high speed HP217 computer to 
control and monitor the echosounder. This enables us to 
obtain real-time data collection and analysis with the 
ability to store and reproduce the atmospheric profile 
probs at will. This information should help quantify the 
effects of low level turbulence on electro-optical systems 
performance. lt may then be possible to negate these 
turbulent eftects by Simply mounting sensors) and 


transmitters well above the turbulent layers or by adaptive 


Optics. 
times tnesis discusses the development of a high 
frequency, 5KHz, HeoOus ile echosounder, the criteria 


involved in the actual designs of both the acoustic array 
and the acoustically absorbent enclosure, and the selection 
and function of all the associated component parts of the 
echosounder system. Finally the operational success of the 


system will be discussed and the output data analyzed. 


II. BACKGROUND 


Acoustic echosounders probe the atmosphere by 
transmitting a pulse of acoustic power which is 
subsequently scattered back from the atmosphere by 
temperature and velocity inhomogeneities. The echosounder 
(echosonde) equation, often referred to as the radar 
equation in meteorology, is used Lo determine the 
backscattered acoustic power. This equation is summarized 


by Neff in Reference 10 and is based upon the work of 


Tatarski [Ref. 11] and Little [Ref. 5]. 
Pre = Er Prersji exp(-2ak) Fi com@Rk,f) 9 seul tia Greme 


where 

PR is the electrical power returned from a range R, 

Pr is the electrical power transmitted at frequency f, 

RR is the efficiency of conversion from acoustic power 
to electrical power by the transducer, 

Et is the efficiency of conversion from electrical 
power to acoustic power by the transducer, 

exp(-2aR) is the round trip power loss due to 
attenuation where a is the average attenuation (meters! ) 
to the scattering volume at the range R (meters), 

oo(R,f) is the scattering cross section per unit volume 


at a distance R and frequency f, 
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c is the local speed of sound (meters/second), 

tT is the pulse length (seconds), 

A is the aperture area of the antenna (meters? ), 

Ris the range (meters), and 

Gis the effective-aperture factor of the antenna. 

Kmpirically measuring or calculating the values for all 
other terms, one can use this equation to determine 
oo(R,f), the scattering cross section per unit volume; that 
is, the fraction of incident power backscattered per unit 
distance into a unit solid angle at a frequency f. Based 
upon experimental results, Tatarski [{[Ref. 11] expresses the 
acoustic backscatter cross section per unit volume, 


morn. fh DY, 


Cre 
Pon t= OmOO39 Ach/3. ------- ; 
TOs 
where 
ee = ZY X is the incident acoustic wavenumber at 


wavelength A, 
To is the local mean temperature in degrees Kelvin, and 
Ct? is the temperature structure parameter. 

Combining this equation with the echosonde equation, one 


obtains a volume-averaged measure of CT*?. 


i i! 2 it PR 
CT? = sccecerr eree-- To? k-1/3 eece ---- ---- RF exp( Zak) 
0.0039 ERET er AG Py. 


vi 


Hall and Wescott {[{Ref. 12} calculated a beam-shape 
compensation factor of 0.40 for a piston source antenna 
with a uniformly illuminated square aperture. This value 
is the same as the effective-aperture factor, G@ and can be 
substituted into the above equation. Approximating the 
aperture area of the antenna to be equivalent to 25 times 
the aperture area of a single speaker having a diameter of 
7.620 centimeters, we get a value of 0.1140 square meters. 
Combining these values with the numerical constants in the 


above equation, we can simplify the equation for Cr?. 


iH di PR 
Cr? = 11245 ------ To? k-1/3 ee-=- =--=-=- R2 exp(ZaR) 
ER ET CT Py 
The efficiency of conversion factors for oun 


echosounder were determined with the use of an IVIE 
Hlectronics IE-30A Audio Analyzer Sound Level Meter, a 
HP3S3SI4A “Punmetaon Generator, a HP3561A Dynamic Signal 
Analyzer, a reference speaker, the acoustic array and the 
Naval Postgraduate School anechoic chamber. The function 
generator was used to transmit a 1.0 volt, 5.0 KHz, 
Sinusoidal signal through the reference speaker. The 
acoustic array was employed as a receiver and was aligned 
along the centroid of the main lobe of the reference 
speaker signal. With the acoustic array in position, the 


sound level meter was used to determine the total power 


bee 


available for reception by the acoustic array; that is, the 
total power actually ensonified within the solid angle of 
the acoustic array. Ten independent measurements were made 
with the sound level meter and then averaged to yield a 
power value of 73.7 dB. Using a conversion factor and the 
array aperture area determined above, the total acoustical 
power, Pa, available at the acoustic array is given by the 


equation below. 


Pat Of Si eee AOL 27 Watts / més. 0.11401 m2 ) 
eis, 


Pa = 2.67 * 10-S Watts 


The dynamic signal analyzer was then used to measure the 
actual power received by the acoustic array and converted 
into electrical power. Four independent measurements were 
recorded and averaged to yield a value of 4.035 mVrms. The 
acoustic array impedance was previously measured and 
determined to be 12.1 +0.4 Q. These values were used in 
the following equation to determine Pe, the total 


electrical power converted by the acoustic array. 


(ajo wee O= >) Vrms) 
Pe = seer ee eee e senor seer e----- = 1.34 © 10-6 Watts 


The efficiency of conversion factor for acoustic to 
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electrical power is the ratio of Fa to Fe ywor, 


The value for the efficiency of conversion for electrical 


to acoustic power, ET, is not readily calculatable with the 


equipment at hand. However, based upon the design of the 
speakers used, the values for ER and Et Should be 
approximately equal. Furthermore, allowing if Ox 


experimental error especially in the reading of the sound 


level meter, the best accuracy for the values of ER and ET 
is 0.5 +0.1. Substituting conversion factor values of 0.5 
into the above equation, we can further simplify the 


equation for CT*. 


Ros ER 
CT® = 44980 k-1/3 -w2-2-- ---= RF exp(2aR) 
CT PT 


The above equation was incorporated into a computer 
program that provides the system control, data acquisition 
and data processing techniques [Ref. 3]. The reduced data 
was then used to measure the temperature structure 
parameter as a function of time of day and altitude for 


various sites. 
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III. SYSTEM DESIGN AND EQUIPMENT DEVELOPMENT 


A. SPEAKERS 

In the design of our echosounder, it was determined 
that a rapid decay time was required to obtain accurate 
short range information. Weight restrictions involved with 
equipment transportation plus high efficiency in the 2KHz 
to SKHz range led to our decision to use piezo ceramic 
speakers. The Motorola KSN 1005A speaker was selected 
based on these requirements and the specifications charted 
in the Motorola catalog [{Ref. 13} and reproduced below in 


Figures 1 through 3. 


Nominal Power, Impedance, and Distortion Ratings 
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“Fig. 1. Speaker Ratings 
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Fig. 3. Speaker Dimensions 
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These charts indicate that a maximum response for our 
speakers occurs at a resonant frequency of 5000 Hertz. 
This frequency was used as the baseline from which all our 
measurements were made. 

Using the speakers in an anechoic chamber, the average 
e-1 voltage decay time was measured to be approximately 900 
usec (Fig. 4). This decay time translates into a sound 
propagation distance of just over 15.0 centimeters (at STP) 
from the speakers. Considering our requirements, this 


speaker was ideally suited to serve our purpose. 





Fig. 4. Speaker Decay Time Trace 
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B. ACOUSTIC ARRAY 

Our next consideration was the echosounder array 
pattern. Ideally, the acoustic sources should be placed 
exactly one half wavelength apart. At a frequency of 5000 
Hertz, this would require spacings of 3.4 centimeters (at 
STP) which is physically impossible for the speakers we 
have chosen. The closest possible spacing was 7.62 
centimeters between sources after shaving off the flange of 
the horn (Higa = 


From Kinsler, et al. [Ref. 14: p. 195], the equation 





for the directionality factor of a simple line array is 


derived as: 


N 
sin ( --- kd sin 6 ) 
1 2 
H(N, 9S) Bm pce a eee een: ; 
N 1 
sin ( --- kd sin ® ) 
a 
where 


k is the wave number (2n/\X ), 

dis the distance between sources, 

N is the number of sources, and 

8 is the angle measured from a line perpendicular to 
the array to the direction of interest. 

However, this equation assumes simple point sources 
which does not accurately describe the speakers we have 


chosen. It was necessary to couple this equation to the 
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directionality factor for a piston source which is also 


Identified in Kinsler, et al. [Ref. 14: p. 108]. 


2Jd1 ( ka sin 8 ) 


where 

k 1s the wave number, 

ais the radius of the piston source, 

@ is the angle measured from a line perpendicular to 
the array to the direction of interest, and 

Ji is a first order Bessel function. 

These two equations were combined to produce the 
equation of directionality for a linear array of piston 


sources, L(N,8), by simple multiplication. 


hemo) = HON, Gy Dice) 


This equation was incorporated into a computer program 
[Appendix A} and used to generate linear array patterns for 


varying numbers of piston sources (Figs. 5 through 9). 
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Fig. 5. Three Element Array 


-10 


DECIBELS 


PATTERN FOR 






PATTERN FOR 4 ELEMENTS 


-18 


DECIBELS 
ru 
G 


! 
GJ 
q 





Fig. 6. Four Element Array 
5 ELEMENTS 






Fig. 7. Five Element Array 
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Fig. 9. Seven Element Array 


It was concluded from these plots of the linear array beam 
patterns for varying numbers of piston sources that five 
linear elements would provide the best combination of 
LOrverd = aereclLionality, Sidelobe suppression, ee cal 
size, and a relatively low cost. Then, in order to enhance 
both array efficiency and symmetry, we settled on a five by 
five element array design with 7.62 centimeter spacing 
between speakers in both the vertical and horizontal 
directions. 

After we verified the manufacturer’s polarity 
designation for 35 speakers and inspected each of these 
speakers for physical defects and loose wiring connections, 
we obtained Lissajous plots for each individual speaker. 
Figure 10 represents a typical Lissajous plot for one of 


the speakers incorporated into the five by five element 


array. Based upon the speaker’s output to input voltage 
ratios as illustrated by these Lissajous plots, we were 
able to rank all of our speakers by signal efficiency. te 


was based upon this criteria that we selected the 25 most 
efficient speakers for use in the five by five element 
array, placing the best speakers at the center and 
subsequent ranking speakers further toward the sides and 


G@rrers . 
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Fig. 10. Lissajous Plot 


The 25 selected speakers were mounted in a five by five 
planar array on a balsa wood insulated bilayered sheet 
metal board. After wiring all the speakers in series, we 
surrounded all the electrical connections and speaker backs 
with two 3.0 centimeter layers OF foam insulation 
sandwiching a 1.0 millimeter lead sheet. Then the entire 
array mounting was enclosed ina 44 by 44 by 5 centimeter 
Sheet metal box. This design was chosen to suppress 
virtually all acoustic energy propagating out the rear 
hemisphere of the array, while shielding the array from any 


external electrical interference (Figs. 11 and 12). 
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technique 


useful 


atmosphere. [Ree COndemE omit ltze | the acoustic waves 
transmitted and later received by this remote sensing 
method efficiently, it is essential to have an efficient 
antenna with highly directive beams and strongly iponessed 
Sidelobes. Antenna design becomes increasingly more 
important in a noisy environment where noise pollution 
within the sidelobes may dominate the desired signal within 
the main lobe. Hall and Wescott [Ref. 12] showed that 
sidelobe suppression improved with higher frequencies. 
Their studies showed that the measured 90 degree sidelobe 
suppression ranged from 38 dB at 1 KHz to 50 dB at £ “5 KHz. 
Furthermore, any Significant improvement in sidelobe 
Suppression could only be obtained by surrounding the 
antenna with an acoustic energy absorbing cuff or shroud. 
In an effort to maximize our antenna main lobe to 
sidelobe power ratio, we intend to operate at  5KHz. 
Additionally, we have designed an acoustic energy absorbing 
enclosure. Many designs were considered based upon 
previous research in the field of echosounding [Refs. 15 
through 18}. In addition, we obtained the actual acoustic 
beam patterns for our array using a computer program 


written by LCDR Butler [Ref. 19] which we modified for our 


purposes. This modified version of LCDR Butler’s computer 
program is in Appendix B. By rotating the array in an 
anechoic chamber, we were able to produce highly accurate 


polar plots of the array beam patterns (Figs. 13 and 14). 
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Fig. 13. Polar Piot of 5 X 5 Array Beam Pattern 


26 


DIAGONAL 5 X 5 ACOUSTIC ARRAY BEAM PATTERN 
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Polar Plot of 5 X 5 Array Beam Pattern 
at a 459 Aspect 
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A comparison of the theoretical beam pattern for the 
five element, piston source array in Figure 7 with the 
experimentally measured polar plot of the 5 X 5 acoustic 
array in Figure 13 and the associated coordinates 
corresponding to the individual data points illustrate 
Significant agreement. As predicted by our computer model 
{Appendix A], the actual acoustic beam pattern is comprised 
of nine distinct lobes in the forward hemisphere. Both the 
positions of these lobes and their relative amplitudes are 
Within 2.0 percent of their predicted values. The computer 
modeled main lobe full width angle of 19.5 degrees very 
closely predicted the experimentally measured main lobe 
full width angle of 19.0 degrees. Reflected, scattered, 
diffracted and electronic noise all contribute to the 
clutter near the origin of Figure 13. This is the primary 
reason why the minima don’t actually go to zero as 
predicted by our model. Additional reasons include 
mechanical noise in the device used to rotate the acoustic 
array and vibrational noise in the array itself. The 
apparent lobe in the upper left hand quadrant is an anomaly 
seen in all polar plots herein, and corresponds to acoustic 
energy reflected off the door and window of the anechoic 
chamber. Figure 14 represents the acoustic array mounted 
at a 45° aspect. A computer model for this orientation is 
a very complicated task and was not performed at this time. 


However, this polar plot of the acoustic array beam 


Zo 


Pettern, and in particular, the appreciable sidelobe 
suppression strongly support the need for further array 
design studies. 

Since it was our aim to suppress all FP dailepes and 
utilize the main lobe, we chose not to taper our enclosure 
as most previous researchers had. Rather we designed the 


enclosure based upon the dimensions of the array itself and 


the acoustic beams it generated. The height of the 
enclosure is 160 cm. Allowing for the base foam and lead 
insulation, the array housing and the base framework 


itself, the actual array piston sources are positioned 143 
em from the mouth of the enclosure. The length and width 
of the enclosure are both 65.) cm. Again allowing for 
construction materials, the inner dimensions are both 48 
cm. Hence, the position of the center speaker is 24 cm 
from any of the four sides. These dimensions produce the 
unobstructed 19.0 degree path into which the main lobe is 


transmitted. 


fea 24 em. 7 t4aeem ) = 19.050 = 19.09 


BPivweod was Used for the construction of the enclosure 
and provided not only a rigid, inexpensive framework, but 
also provided an impedance mismatch which helped to 
attenuate external noise interference. Anticipating all 
kinds of weather conditions during data collection, the 


plywood enclosure was first waterproofed with four coats Gir 


marine varnish. Grooved joints, caulking and weather 
stripping were also design considerations. 

A millimeter layer of lead can suppress an acoustic 
Signal by reflection as much as 40 dB (Figs. 15 and 16). 
About 7.0 centimeters of convoluted foam can suppress a 
Signal another 3 to 4 dB (Figs. 15 and 17). The primary 
purpose of the convoluted foam is to attenuate the 
sidelobes which repeatedly reflect off the enclosure 
interior at glancing incident angles. Together the lead and 
convoluted foam make an extremely efficient absorbing 


material for use in our enclosure. 
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Fig. 15. Signal Response Reference, No Insulation 
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Fig. 17. Signal Suppression by Foam 
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Fig. 19. Enclosure Photograph, Interior 
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D. ADDITIONAL HARDWARE 

lig) -Sielolis re Weyl to the acoustic array and enclosure 
discussed in the previous sections, there are primarily six 
more pieces of equipment essential to the func es the 
echosounder. Below is a brief description of these 
additional components. 

1. HP200 Series Computer 

The Hewlett-Packard (HP) 200 Series Computer 

includes a 20 megabyte hard disk drive, a floppy disk drive 
and an associated printer and monitor. The HP2Z00 Series 
Computer is the central control component for the entire 
echosounder arrangement. The computer used is an HP21i7 
programmed in Basic 3.0 and equipped with an Infotek BC203 
Basic Compiler and an inteotrexs HEZI0™” Kioating Point 
Accelerator to enhance the speed of program execution. The 
computer is used to maintain timing and process control, to 
define all parameters for the HP3314A Function Generator as 
well as to execute the trigger command that produces the 
echosounder transmitted acoustic pulse. The computer also 
processes the data from the acoustic array via the pre- 
amplifier, bandpass filter and analog to digital converter. 
The computer then conducts data reduction routines to 
produce and display a high resolution atmospheric profile. 

2 meee IAA Nunction Generator 

The Hewlett-Packard (HP) 3314A Function Generator 


is a multimode function generator capable of providing 


a0 


Sine, square and triangular wave functions as well as any 
desired waveform ranging in frequency from 0.001 Hertz to 
19.999 Megahertz. The HP3314A Function Generator was used 
to supply a pulse of an integer number of sinusoidal cycles 
of constant amplitude to the QSC Model 1700 Audio 
Amplifier. Typically a pulse of 100 sinusoidal cycles was 
used. A constant frequency setting of 5000 Hertz was used 
for all data runs. 

3. @5C Model 1700 Audio Amplifier 

The QSC Model 1700 Audio Amplifier is a high power 
amplifier which can supply 350 watts over an 8 ohm load. 
This amplifier was used to boost the output voltage of the 
function generator by a factor of 20 from 1.5 volts to 30 
volts: The power supplied to each speaker in the array 
during operation was then 37.1 watts as indicated in the 
equation below. 
( Voltagerms )@ ( 30 )2 
EOS a ee = 37.1 Watts 
Impedance 2 ee lc oan 
4. Pre-amplifier 
The pre-amplifier was designed and constructed by 

Walters [Ref. 20] to supply a gain of 1000 to the returned 
Signal and isolate the data aquisition system from the 
transmitted pulse. The pre-amplifier may be thought of as 
a safety and switching mechanism for the system. An om 


1007 and an LT 1037 Operational Amplifier were selected for 
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use in the pre-amplifier based upon their low noise 
properties which were evaluated by physically incorporating 
them into the device and measuring their noise and gain 
characteristics with a HP3561A Dynamic Signal Analyzer. 
The IN 4000 Series Rectifiers isolate the power amplifier 
from the receiver between pulses and limit the voltage 
applied to the operational amplifiers. The electrical 
diagram of this device is outlined in Figure 21. 
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Fig. 21. Schematic of Pre-amplifier 
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5. Rockland Wavetek Model 852 Bandpass Filter 


The Rockland Wavetek Model 852 Filter operates as a 
48 dB per octave bandpass filter that suppresses broadband 
noise in the system. High and low bandpass settings of 
5500 Hertz and 4500 Hertz respectively were used with a 
roll off of 48 dB per octave. The filter response at these 


settings is illustrated in Figure 22. 
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Fig. 22. Bandpass Filter Response 


The 3dB (half power) bandwidth was approximately 1 KHz 
wide. Although it would have been desirable to reduce the 
bandpass interval to between 50 and 100 Hertz without 


impairing signal response, the Rockland Wavetek Model 852 
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Bandpass Filter was the best filter available. ie 
performed adequately despite the large 1 KHz bandpass. 
6. Infotek AD200 Analog to Digital Converter 
An Infotek AD200 12 bit Analog to Digital Converter 
was used to sample the amplified signal voltage and 
digitize it for the computer system. A typical sampling 


frequency used was 20000 Hertz. 


The arrangement of the above equipment is illustrated — 


mmr ieure 23. 
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Fig. 23. Echosounder Equipment Arrangement 
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BH. SOFTWARE 

The software and the HP217 Computer were responsible 
for controlling and monitoring every operational aspect of 
the system. LT Wroblewski [Ref. 3] provides an in-depth 
discussion of the echosounder software which is also 


summarized in Appendix C. 
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IV. DATA ANALYSIS 


A. ECHOSOUNDER PERFORMANCE 

An analysis of the acoustic echosounder output was 
conducted to determine the validity of previously 
determined echosounder parameters. The e-1l1 decay time 
constant for the applied voltage of a single speaker was 
previously calculated to be 900 us which corresponds to a 
distance of 15.3 cm. The e-1 decay time constant for the 
applied voltage of the entire acoustic array was determined 
by the same experimental procedure to be 1060 us which 
corresponds to a distance of 18.0 cm. Although the 25 
speakers in the acoustic array transmitting in unison 
should have approximately the same e-1 decay time constant 
as the single speaker, it is not surprising that the array 
e-! decay time constant is larger. Using ten time 
intervals, an expectation of 99.99 percent decay of the 
transmitted signal should occur within a distance of 180.0 
cm or just under 2.0 meters. Reviewing typical echosounder 
traces and Ct2 plots {Appendix D], it is evident that such 
BeeeiIOnimrecovery time 1S actually quite accurate. The 
recovery distance is typically found to be on the order of 
oometers or 53 ms. This “blind zone" is not uncommon in 
most radar systems and is in part a return signal from 


ground clutter surrounding the acoustic array and its 
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enclosure immediately following pulse transmission. An 
analog to digital intensity factor of 4000 is assigned to 
the heavily shaded region of the recovery zone. Since the 
computer software has generated this value as the product 
of the actual analog to digital “count” value and the range 
from the acoustic echosounder, this value of 4000 must be 
divided by the range at the recovery zone boundary to give 
us an actual analog to digital “count value. This result 
is then multiplied by an experimentally measured conversion 
factor of 3.34 ° 10-8 voltsrms/count. This result may then 
be used in the determination of the total electrical power 
received and converted by the acoustic array. The value 
for the total electrical power received by the acoustic 


array is calculated in the equations below. 


Voltagerms = (3.34 © 10-8v/count) (4000mecounts) (Y9m)-1 


Voltagerms = 1.48 * 10-5 volts 


Power = (Voltagerms )2 ¢ (Array Impedance)-1! 


(1.48 + 10-8 volts)? + (12.1 Q)-1 


iI 


Power 


Powerw = 1.82 = 10°21 Watts 


Using this received electrical power value and the known 


transmitted acoustic power value of 37.1 Watts, we can 
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solve for the e-1! decay time constant of the applied 


voltage. The product of the total acoustic power received 
and the efficiency of conversion factor ( ER = 0.50 ), 
yields the total electrical power received. Hence, the 


total acoustic power received is 3.64 - 10-11 Watts. 
Although the distance of 9 meters corresponds to a time 
interval of 53ms, this value is actually inflated by the 
system pulse length of 20 ms. Compensating for this, the 
actual recovery zone time interval is 33 ms. The e-2 decay 
time constant of the power is equal to the e-! decay time 
constant of the applied voltage. These values are 


calculated below. 


Tm www ww enw ww wee ewe ewe =e, 20 mS ’ 
tm (6.37 © 10-12) 


This e-1 decay time constant of the applied voltage 
corresponds to a distance of 43.5 cm and is about twice the 
previously determined value. 

Consideration was also given to the construction of a 


hexagonal planar array rather than a square planar array 
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based upon the astounding sidelobe Suppression for the five 
by five element, square array at a 459 orientation as 
illustrated in Figure 14. To investigate the potential of 
such a close packed design, two small seven element 
hexagonal arrays were constructed. The first array design 
utilized a close packed spacing of 7.62 cm, while the 
second array design employed a 3/2 wavelength or 10.2 cm 
spacing. Polar plots of the beam patterns for both these 
hexagonal array designs were made by rotating the arrays in 
an anechoic chamber with the main diagonal both horizontal 
and vertical. In both cases the first sidelobe suppression 
for both close packed orientations was down by nearly 15 
dB, while the first sidelobe suppression for both 3/2 
wavelength spacing orientations was only down by about 7 
abe Furthermore, the ratios of the size of the main lobe 
to those of the sidelobes is much greater for the close 
packed design than they are for the 3/2 wavelength spacing 
arrangement for bDoun the vertical and horizontal 
orientations. The polar plots of the array beam patterns 
for both close packed and 3/2 wavelength spacing hexagonal 
designs are illustrated in Figures 24 through 27. These 
figures strongly warrant further study of the hexagonal 


close packed array design. 
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HORIZONTAL DIAGONAL, CLOSE PACKED, 7 SPEAKER HEXAGONAL ARRAY 





R=3. 60 M MIKE VT= 5.6364 INPUT=5.0 V FREQ=5. 0 KHZ 


Fig. 24. Polar Plot of Close Packed Hexagonal Array 
(Horizontal Diagonal) 
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VERTICAL DIAGONAL, CLOSE PACKED, 7 SPEAKER HEXAGONAL ARRAY 





R=3.60 M MIKE VI= 5.1559 INPUT=5.0 V FREQ=5. 0 KHZ 


Fig. 25. Polar Plot of Close Packed Hexagonal Array 
(Vertical Diagonal ) 
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HORIZONTAL DIAGONAL, 3/2 WAVELENGTH SPACING. 7 SPEAKER HEXAGONAL ARRAY 





R=3.60 M MIKE VT= 5.18 INPUT=5.0 V FREQ=5.0 KHZ 


Fig. 26. Polar Plot of 3/2 Wavelength Hexagonal Array 
(Horizontal Diagonal) 
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VERTICAL DIAGONAL, 3/2 WAVELENGTH SPACING, 7 SPEAKER HEXAGONAL ARRAY 





R=3.60 M MIKE VTI= 6.0215 INPUT=5.0 V FREQ=5. 0 KHZ 


Fig. 27. Polar Plot of 3/2 Wavelength Hexagonal Array 
(Vertical Diagonal) 
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B. SITE EVALUATION 

Echosounder data was collected at two different 
locations. The primary data collection site was the upper 
roof of Spanagel Hall at the Naval Postgraduate Sehere 
Monterey, California. This site was chosen simply for 
convenience of data collection. Data gathered at this 
location is believed to represent the California coast 
during the spring near sea level. The second data 
collection site chosen was in the vicinity of the 24 inch 
telescope at Lick Observatory, San Jose, California. This 
Site is located atop Mt. Hamilton at an altitude of 
approximately 5700 feet and nearly 20 miles inland from the 
coast. 

These two data collection sites represent areas of 
differing atmospheric air pressures, water vapor pressures, 
local temperature ranges, and local wind velocity ranges. 
These characteristics akl play important roles in effecting 
the local atmospheric turbulent conditions and thereby the 
atmospheric structure parameter, CT2. 

In addition to collecting echosounder data at hick 
Observatory, simultaneous measurements of the isoplanatic 
angle (80) and spatial coherence length (ro) were made with 
systems developed by Walters [Refs. i ance |e A basic 
knowledge of these two systems is necessary to understand 


the correlation procedures made. The isoplanatic angle 
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(Go) 1s primarily an upper atmospheric measurement which 
indicates atmospheric disruptions at a range of two to 15 
kilometers above the Earth’s surface. On the otherhand, 
the spatial coherence length (ro) is auniform measure of 
the effects of the entire atmospheric blanket on coherent 
light transmission. A cloSe comparison of the acoustic 
echosounder data, the isoplanatic angle (80) measurements 
and the spatial coherence length (ro) measurements should 
give us an indication of the turbulent contributions from 
the lower and upper troposphere as well as the stratosphere 
above. 

Based upon isoplanatic angle (80) and spatial coherence 
length (ro) measurements at Mt. Wilson near Los Angeles, 
California, a strong correlation between the isoplanatic 
angle (80) and spatial coherence length (ro) measurements 
occurs if the low altitude boundary layer contribution is 
sufficiently small. This strong correlation helps to 
reinforce the overall description of the atmosphere at the 
time of data collection. Figures 28 and 29 graphically 
illustrate the atmospheric measurements made at Mt. Wilson 
on 2 April 1987 by Walters. The strong correlation between 
the isoplanatic angle (80) and the spatial coherence length 
(ro) is especially evident between the hours of 0/700 and 
1300 universal time. The close tracking of these two 


measurements during this time interval indicate that the 
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Fig. 28. Isoplanatic Angle Measurements, Mt. Wilson 
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Fig. 29. Spatial Coherence Length Measurements, Mt. Wilson 
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upper atmospheric eenduLrons, as measured by the 
isoplanatic angle Bon). are dominating the entire 
atmospheric profile as measured by the spatial coherence 
length (ro). Unfortunately, the existence or non-existence 
of any turbulent surface effects cannot be ascertained by 
the employment of the above two systems alone. However, 
employment of these two systems together with the acoustic 
echosounder should enable us to produce a complete 
atmospheric profile with strong correlation between all 
three atmospheric measurements. 

On 9 and 10 April 1987, all three systems were operated 
at Lick “ObpSernrvacor,- Again good correlation between the 
isoplanatic angle (80) and the spatial coherence length 
(rom measurements was noted. Im addition strong 


correlation between the spatial coherence length (ro) 


measurements and the acoustic echosounder data was 
conclusively present. A comparison of the isoplanatic 
angle (Go ) and the spatial coherence length (ro) 


measurements in Figures 30 and 31 shows good correlation 
between the two atmospheric parameters especially during 
the 09307 to 1130 time interval on 9 Aprem 1337 However, 
during subsequent hours the isoplanatic angle (Ho) 
measurements remain high (approximately 12 wurad) indicating 
relatively calm turbulent conditions in the upper 


atmosphere, while the spatial coherence length (ro) values 


se 


drop sharply after 1100 universal time tAGheavi ne dominant. 
and increasing lower atmospheric turbulence. This increase 
in the low level turbulence should be evident in the 
echosounder data commencing around 1200 universal time 
(0400 local standard time) on 9 April 1987. A comparison 
of the echosounder data in Figures Sa stosoeugh 34 


illustrates this increase in the local surface turbulence. 
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Fig. 30. Isoplanatic Angle Measurements, Lick Observatory 
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Fig. 31. Spatial Coherence Length Measurements, 
Lick Observatory 
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Fig. 32. Echosounder Trace, Lick Observatory 


o4 


Pie KOS SeRy Fil OR 


138.2 


=) alelt- 


fos 


“ss Neo rane - 
eee “ee ele + pet mee ote so oe 6 © ser Nsegpars N*. - 
Pye aggete MBSE. AN RGD ol Se ERE Ra et} 
Beet ee: eas 0) SA aineers “eo ote iets ogee Mee 8 
: 


4:15 






cece we es Bees VAN AE GL eB ‘ 
: woe oML, ONS ets tk its anette: ee 
Caran Seite 238 Bl ae ss EE Salama co 
> oes a eeens 22*.20° : : : eee *. ° ets Smet sites “eRe » ‘ 


° 
eee ete ee me 


an 
=enee 


peers Abie: 
MLL Fees, Sie se 


ee BMC ES on. 


eo ge RN at: ee a ar SU See 
Br ecto eaca fae eS ofed Ee ae a 
ewedecuee seca eerie Sica eogeren sur ©. sen PSs ee oa 
oo ee 
° ° ee eee ° ae 
ged? xr 
xine 


ce eet eeeee 


eet TS 
. 


¢ 
4 


ese @8 






Sd3L3aW NI JONNY 


MGS. Leis Mesa ep) 


Lick Observatory 


Echosounder Trace, 


ae 


Fig. 


Piero noer VA TOR y 


13 


° 
oe e208 


ese © @ © @ 
. 


Q 
= 


SYSLSIW NI 3oruYd 





1@ 


Ww 


( Sae)) 


cre 


ieee 


Lick Observatory 


Echosounder Trace, 


34. 


Fig. 


30 


The strong correlation between the echosounder data around 
0400 to 0500 local standard time and the spatial coherence 
length (ro) measurements at 1200 universal time, combined 
with the lack of correlation between the isoplanatic angle 
(90) and the spatial coherence length (ro) measurements 
indicate that the lower atmospheric and surface turbulent 
layers are dominating the atmospheric profile during this 
time period. 


Data collected on 10 April 1987, again illustrate the 


strong correlation between the three atmospheric 
measurements made. A comparison of the isoplanatic angle 
(S80) and the spatial coherence length (ro) measurements 


during the time interval of O700 and 1300 universal time 
indicate steady turbulent conditions in the upper 
atmosphere and greatly varying turbulent conditions at 
lower atmospheric levels. This is evident in Figures 35 
and 36 by the consistent values of the isoplanatic angle 
(Qo) during the time period compared with the steady 
increase and eventual decline of the spatial coherence 
length (ro) values during the same time interval. The 
trace variations in Figure 36 during the hours of 0800 and 
1200 universal time are indicative of a period oat 
diminishing lower atmospheric Or surface turbulence 
followed by the onset of an increasingly turbulent period 
around 1200. This turbulent trend MS. “Gils ime ht Ve Ly 


supported by the echosounder data in Figures 37 through 44. 
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The large 200 mm spatial coherence lengths around 1100 
universal time are consistent with the low turbulence 
evident in the echosounder profiles around 0300 local 
standard time. 

Only acoustic echosounder measurements were made on the 
upper roof ot Spanagel Hall at the Naval Postgraduate 
Schou Data runs on 26 and Zi April yigee include both 
echosounder data and the associated Bér2 plots. Ox 
particular interest ais the plovpmer eee mon the 2Z/th of 
April in Figure 45, which shows the maritime boundary 
inversion layer around 100 m being perturbed by convective 
plumes at lower altitudes. ‘Sixteen additional echosounder 
output traces were included in Appendix D to exhibit 
typical Atnegpiewmwemmpaetivity. Many of the plots, such as 
the 14:45, 153300 and the 17:30 of ZemAeri1! and thew ie 
Li:15,. 23:00, 242255 Vat Gs icin Cle ty tomy eee 
April, have convective plumes which are prevalent whenever 
a heat flux between the surface and atmosphere exists. The 
plots of 18:15 and 19:15 on the 2Zetn of e8ent!) and Ct wee ie 
18:45, and 20°00 on the @27/th of Apmi welea weashew the 
passage of the neutral event which is encountered when the 
atmospheric amd  surtace Gempem@turcweet se neice DeCOMme—s 
negligible. “Finally, tne plot Of 6 750 er 2 ompinee) Can oe 
associated with strong winds which exhibit a somewhat 
iniform return for all altitudes across the entire 15 


minute interval of the trace. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A high frequency, 5kHz, phased array echosounder 
constructed to measure low level turbulence appears to work 
within the 10 to 135 meter altitude range. When using the 
device, detailed profiles of the short range atmospheric 
density fluctuations were obtained. The protiles were 
found to correlate very well with the measurements of the 
isoplanatic angle (80) and the spatial coherence length 
(ro) during periods of simultaneous operation. This “shone 
range echosounder, when used in conjunction with the other 
atmospheric measuring devices, is an invaluable tool. The 
graphical Oultpue provides us with a more complete 
description of the atmosphere that can be used to calculate 
ni the atmospheric index of refraction structure 
parameter. 

Further software developments such as the incorporation 
of Doppler and Fast Fourier Transform routines into the 
echosounder program will determine the radial velocity 
profile of the return signal. Such an improvement will 
determine the velocity of the probed air masses’ and plot 
these echosounder traces as a function of color intensity. 

Other areas of further reasearch include the design and 
testing of different array patterns. One such pattern, the 


hexagonal array, is already undergoing tests. As the 


64 


Speaker array evolves, it is also evident that the 
enclosure must follow suit to accomodate the new array 
design and associated beam patterns. Further improvements 
in software routines are also inevitable. It is always 
desirable to store on a disc all data collected at a site. 
Presently, the floppy disc capacity is inadequate for the 
storage of more than two hours of data. Furthermore, if 
FFT and Doppler routines are to be added, the computational 
speed of the computer will inhibit the pulse repetition 
rate. A special data acquisition technique called "Direct 
Memory Access’ may then have to be introduced into the HP 
computer to allow Simultaneous data Goiuleetion. and 


processing. 


Finally, thas Produed may be used in various 
- applications not already addressed. One such use may be to 
measure windshear at airports. Other researchers have 


expressed interest in using the echosounder to measure 


ane tC atmospheric Goneitions auring meteorological 
surveys. Assuredly, as the device evolves into a more 
compact and highly efficient instrument, its range ot 


Sepiication will continue to expand. 
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APPENDIX A 
LINEAR ARRAY PROGRAM 


{ANTENNA BEAM PATTERN FOR AN N ELEMENT LINEAR ARRAY 

GCE Ar 

Ri=.Q3581 

Fr=5200@ 

Ka=(Ri*2*Pl*Fr )/343 

INPUT “ENTER THE # OF DRIVERS ~ N 

D=.@762 

Kd=(2*PI*D*Fr )/343 

Kd2=Kd/2. 

Num_decades=4 

Scale=1@*Num_decades 

Thetastep=.2@@1 

GCLEAR 

VIEWPORT 20,820 ,22,3@ 

WINDOW =2051905 =2-17¢ | 

AXES 18@8/(2*Num_decades ) ,1@@/(2.*Num_decades) ,@ ,@,28,2@ 

MOVE 128,103 

LAGE FAn en teFOR “SN (BSE n toe 

CS ire 

FOR = \aees 

MOVE 1S C2 See) | 
M=C1O=( (B=) 41a )) 
LABEL M 

NE Kelas 

LOTR eee 2 

MOVE — |S e40 

LABEL DEGIBEES © 

LDIR @ 

FOR Theta=Thetastep TO PI/2 STEP Thetastep 
Arg=Ka*SIN( Theta) 
Aroc=22FNU | AreD 
Arg2=Arg2/Arg 
Kds=Kd2*SIN( Theta) 
Htheta=SIN(Kds#N)/(N*SIN(Kds )) 
Htheta=Htheta*Arg2. 
Htheta=Scale*Htheta*Htheta 
Htheta=LGT(Htheta)/Num_decades 
X=Htheta*SIN( Theta) 
Y=Htheta*COS( Theta) 

PLOT X#*1@@,Y*1@@ 
!PRINT Theta,Htheta 
NEXT Theta 


END 
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DEF FNJ1(X) 

Pemeloset PUG PoOmmicmiar. 19642 DLW 

hem Oe (HEN 60 One to 

X3=X/3. 

X3=X3*X3 

Px=.54X54(-.56249985+xX5*( .21095575+X3*(-.03954289+xX3* 
( ,@Q@443319+X3*(-.@0031761+X3*.000@011@9))))) 

J{=Fx*X 

Rete RNe 


it sm Noe a1 X 


Fx=.79788456+X3*( .00000156+X3*( .@1659667+X3*( .@001710S+X3* 
(—.@@2495114+X3*( .@@1 1 3653-X3*.@0020@33))))) 
Tx=X-2.35619449+X3*( .12499612+X3*( .@@Q@S6S5+X3*(-.00637879+ 
X3*( .Q@0Q074348+X3*( .00079824-X35*.00@29166))))) 

J1=Fx*COS(Tx )/SQR(X) 

RETURN JI 

FNEND 
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APPENDIX B 
LCDR BUTLER’S BEAM PATTERN PROGRAM 


10 PRINT “** BEAM PATTERN COMPUTATION AND PLOT PROGRAM <x" 
QO PRINT [6 66 eo RRR EKER KEKE KREME RKLR ERK ERE ERR KRENEK KE" 
SU 

40 ! 

20 PRENYT “===xx CR <== 

50 ! 

(oa 

30 PRINT CHRGaGaD 

SU PRINT "THIS PROGRAM CONTROLS THE FOLLOWING APPARATUS" 
100 PRINT “TO MEASURE, DISPLAY, PLOT AND STORE THE" 

110 PRINT "DIRECTIONAL CHARACTERISTICS OF A SINGLET," 

120 PRINT “DOUBLET, TRIPLET, AND QUADRUPLET ARRAY.” 


isu 

140 PRINT " INSTRUMENT OE ae 
(el rin lausny i 
i560 PRIN] | Beme47e) HUETINETER Gy cies 724 i 
170 PRINT Hees47 o MUTT IME TER Gieop as ‘ 
a ea " DATA AQGUISTTIUN/CONTROL (03 _ 
130 

200 DISP “TSR OI Ee NS UREN Ge 

omens! | 


220 REM *«** PRESET AND INITIALIZE INSTRUMENTS*** 
230 GUIPUT ROS ee oral 

240 QUIPUT 723; "F3Rsli° ! FREE KUMMSORWniis 

250 OUTPUT 724 ;"F2RATI” ! FREE RUN AUTORANGE AC VOLTS 
260 ! 

2/0 DISP “TNA Zee O Givin. 

26 lame 

290 SU? | LUN ee 

300 CUM ANGLE(150) ,CALC_ANGLE(150) ,NAMES(80] 

(il DE 

S205 7k LN VER 

350 INTE Gales 

340 DIM R¢360) ,VULT(360>) , URNS DATACSGD) Xx C2607 3G oo 
Seay seas 

J6U (MASS SIURGCees See wali 

370 DISP “NUMBER UR DATA POINTS Seep 
JOUREEEP i Uoen) 

3399 INPUT N 

400 DISP “ENTER DESO DG iis: ties is ee 

41Q INPUT NAMES 

420 DISP "WANT NEW DATA CY/N)?" 


68 


430 


780 
730 
300 
810 
S20 
330 
840 
$51 
860 
3/0 
880 
oo) 
300 


Bieter f03 300 

INPUT QS 

IF QS="Y" THEN GOTO 480 

IF Q$="N" THEN GOTO 920 

GOTO 430 

REM *** DEFAULT PARAMETERS **# 


' MUST INSERT FOR LINES 390, 400, 410 AND 810 
f PRIOR TU START UF RUN 

' 

Ccenle NAMES; 360,16 

ASSIGN# 1 TO NAMES 

NAMES="RUN_1" 

FULLSCALE=500 ! LOCK_IN MVAC 

CONV=360/10000 ! CONVERT OHMS TO DEGREES 
See OE TOURS EROURETIEM! PARAMETERS =<=* 
PRINT CHRS (10) 

Peli SeGRearH TITEe: s NAMES 

REM *** MODIFY IN PRECEEDING LINE WHEN NECESSARY *** 
DISP "MANUALLY ENTER DESIRED FREQ FOR SPEAKER DRIVER" 
DISP "ENTER O IF NO CHANGE" 

ise ENTER 1) 70 CHANGE TITLE” 

BEEP 10,600 

Dive ©! CHOICE 

Pees theEN GOTO) 630 

tert lneN GOIO sso 

PRINT CHRS (10) 

ner 709 ;"CLSO0” ' ROTATE Ch 

BpaRTI=1_ 1G N 

Soul y0o-. CLSI ! POWER TO MGTOR UN 

WAIT .5 ! MILLISECONDS 

soe 70S 3; OPNie ? POWER TG MOTOR OFF 

WATT 4000 

miieecco , OHMS DATAKI> 

Ship, 24 - VOLICDD 

poo eels NOW ST 

PRINT “GHMS_DATA=";CHMS_DATA< I> 

ene VElet = {VET <) 


NEXT I 

Pubes) 10 N 

poe, OHMS DATACT), VOLT 1) 
NEXT I 


ASSIGN’ 1 TO »* 

GOSUB PLOT ROUTINE . 
GOTO 910 | 

ENTER NAMES ! NEW TITLES 
GOTO 600 

PRINT “DONE” 


oo 


Sau 
3) 2) 
sou 
340 
950 
360 
a0 
330 
esl 
1000 
1010 
1020 
1030 
1040 


—F arse 
Om 
=O) Ot 
Coo wT 


1080 


mB Ue ee WU = ees ok WA 
NM OWON MAM EWNH— OW 
Coron ooeno°onee ocoem 


END 
DISP "WHAT IS DATA RUN DESIRED?” 
BEEP eso 
INPUT NAMES 
ASSIGN4 1 TO NAMES 
FOR I=1 TON 
READ¥ 1 : OHMS DATACI),VGLTCI) 
NEXT I 
ASSIGN¥ 1 TO = 

GOTO 860 

! 

PLUT ROUTINE: t **=* PLOTS POLAR DATA **~* 
j 


PLOT 1 eismercs 

GRAPH 

GCLEAR 

PEN 1 

CINE aU ies st 
SCALE -60,60,-60,60 
XAXIS 0,10,-40,40 
YAXIS 0,10,-40,40 
CS1 Zev 

LINE Tea 

FOR C= 00040 NS Te ene 
MOVE Oe 
LORG 6 

LABEL L-40 

NEXT L 

LENE See 

CS aaa 

poe 0,30 


! INSERT LABELS 1250 AND 1280 PRIOR TO START OF RUN 


LABEL ;"S X 5 ACOUSTIC ARRAY BEAM PATTERN" 


NCS cee 


MOVE 0,-45 

LABEL ;"R=3.60 M","MIKE VT="\VOLTC1),"INPUT=2.5 ¥"," 
" "PREQ=5.0 KHZ" 

CSIZE 4 

MOVE 45,0 | 

LORG S ! 

LABEL ;"MIKE" 

LINE TYPE 4 

FUR L=10 TO 40 STEP 10 


70 


MOVE ie 0 | | 

mex sl-0 I 360 SEP 5 

Dean =CehS ClDwe=scin Cl) 

NEXT I 

NEXT L 

Steere | 

Sete 

PEN 2 

FOR I=! TO N ! BEGINS BEAM PATTERN COMPUTATION 
IF VOLTCI)>MVOLT THEN MVOQLT=VOLTCI) 
NEXT I 

FOR IT=1 TON 

K=QHMS DATACTI)/10000*360 

K1=K 

Leeds DATACT)<10 THEN K=Kk1+3607295 
IF OHMS_DATACI)>10500 THEN GOTUD 1610 
VECTUR=20*LGT (VOLT¢CI)/MVOQLT)+40 
EINE TYPE 1 

IF VECTOR<1 THEN VECTBR=0 

Pee lek. VECTOR 

X(T) =VECTOR#*COS (CK) 

YCI)=VECTOR*SIN (K) 

eee 1), Cl) 

NEXT I 

Pee 

RETURN 
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APPENDIX C 
ECHOSOUNDER SOFTWARE 


The driving force of this system is found within the 
HP200 Series Computer. The software and computer system 
were responsible for controlling and monitoring every 
operational phase of the hardware components. The program 
which accomplishes this task was entitled "ACRDR" and is 
listed in Reference 3. This program was modeled after a 
program written by Walters [Ref. 20), but each program 
performs a distinctly different computational task. This 
program was written in HP Basic 4.0 and was compiled by an 
Intotek BC203 Basic Compiler to enhance the speed of 
execution. 

The program "“ACRDR" is easily broken down into a number 
of blocks and subroutines which perform specific 
operations. These sections are outlined in a tlowchart 
(Fig. 46). Each block is straightforward in its purpose, 
and the program is designed to be as helpful to the user as 
possible. As a prologue to the actual program code, there 
is a listing of all the program variables with a short 
desription of their use. Such a listing familiarizes a 
user with the computations to be made and also provides 
quick references for any future modifications to be 


implemented. 


ee 


The initialization procedures which set all the 
parameters necessary for data collection follow the 
variable definition section. These parameters control such 
things as the contrast between background and return signal 
and the setting of the computer’s internal clock. All 
options of the HP3314A Function Generator are decided by 
the user at this’ phase. Another initialization involves 
preparing the Infotek AD200 Analog to Digital Converter for 
operation. Finally, internal arrays are dimensioned, 
integer variables are defined and the computer function 


keys are redefined to suit the "ACRDR" program. 


Initialization 

Procedures for: 
HP 3314A 
AD 200 A/D Converter 
Intemal Clock 
Constants 


Zz 
Send C é Plot to the Compute Values of Cy 
Pri e and Plot 
Computation of nnter 


Parameters 
Speed of Sound 
Offset 
Range 
Pulse Length 
Attenuation Coefficient N 


Send Plot to the 
Printer 







Is the 
15 Minute Plot 
Complete? 





Prepare Screen for 
Data Display 


Trigger Pulse from 


HP 3314A Perform Reduction Noise Level 


Collect Return and Compute the Average 


Routines, Store Data ff 
Required 





Fig. 46. Flowchart of Fundamental Program Operations 


T3 


Now that the system 1S ready to operate, certain 
constants must be determined based upon atmospheric 
conditions at the time of data collection. From the 
ambient temperature, a speed of sound in air is determined 


by the equation in Kinsler, et al, [{Ref. 14: p. 106]. 


C= Con eee (ale) en) 


where 

© is the local speed of sound, 

Co is the@epeed of sound in alr at 0 Cele Gsee oo 1 
m/s), and 

T is the local temperature in degrees Celsius. 
Using this value and a user decided sampling frequency and 
data point total, the maximum range of the traces is 


determined by the relation, 


RANGE = C/2 © 1/(SAMPLE FREQUENCY) - N , 


where 

RANGE is the range of the echosounder in meters, 

SAMPLE FREQUENCY is the user desired rate of data 
acquisi Filion wens 

N is the number of data points per trace. 
Typically, the values of many of the parameters are not 


varied. The value ot N above is almost always taken to be 
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the allowed maximum of 16301 data points and the speed of 
sound can be safely estimated to be 340 m/s at room 
temperature. By defining a typical sampling frequency of 
AUOVOSHZ a range of approximately 135 meters is Aerie 

Another parameter determined Ls the atmospheric 
Seeemuation coefficient, tars calculation 1S done in a 
subroutine obtained from Reference 21 and is necessary in 
Zio Dp rorram for the calculation of the temperature 
Structure parameter, CTe. Other computations made prior to 
program execution include the pulse length, wavenumber, and 
the D-C offset of the equipment. 

At this point, the computer is finally ready for data 
eeekeetion. The screen setup displays a distance versus 
time plot and the internal clock of the computer is 
“synchronized with the pulse from the HP3314A Function 
Generator. Following transmission of the pulse, the return 
Signal is received, digitized, stored in an array and data 
reduction commences. A block averaging technique is used 
in which a block of data points the size of the number of 
emitted cycles is summed and averaged. Returns are plotted 
as darkened areas on the echosounder traces with the 
intensity of the darkened area directly related to the 
magnitude of the return signal. This same procedure is 
used in the computation of CT2 with some very important 
ditferences. Aside from Simelyverscorrecting tor nr-é 


spherical divergence, Ct? is corrected for electronic gain. 


ae 


the ratio of power returned to peowempet ranean) ee ene 
efficiency of the speakers as transmitters and receivers, 
the area of the speaker array, the atmospheric attenuation, 
the pulse length, the temperature, the scattering cross 
section per unit volume at a specific range and frequency 
and finally the effective aperture factor of the antenna. 
Additional hyn (Range) is averaged for a particular 
altitude over 15 minute intervals. Finally, after each 
pulse is reduced, a corresponding mean square noise level 
is determined. This is done by averaging each block 
average at maximum range until at least ten values have 
been used in the average. An upper limit of five over the 
average noise level OE habe corresponds to voltage 
fluctuations on the order of U0-' vols) ae) Set pone 
rountine to avoid averaging any strong return signals or 
anomalies such as passing aircraft. After the ten values 
are averaged, every subsequent pulse is averaged into all 
the preceding noise levels and removed from each subsequent 
return signal. 

Ultimately, after each 15 minute interval, or at the 
user’s request, the display terminal image is printed. At 
this time the Cr4¢ computations are conducted and plotted on 
the screen and printer as range versus logio of CT2. Upon 
completion of the printing of these plots. the system 


begins the data processing for the next 15 minute interval. 


16 


There are certain options built into the program to 
allow a user to change various aspects of operation. The 
function keys allow the user to change the local 


temperature, sample frequency or intensity factor during 


program execution. poeitiowallz the wser Gan quit or 
restart the program, print the partial trace on the screen 
or elect to save a trace on a floppy disc. The save 


routine is invoked tor the subsequent 15 minute interval 
atter the appropriate function key is depressed. Saving a 
future trace may be ackward, especially if a user would 
like to keep an interesting trace which is presently on the 
terminal. This problem cannot be readily solved uniess 
each trace is recorded to disc without user intervention. 
At present, this is not done because only eight traces (two 
maeeurso ot Gata) can be written to a floppy disc before it is 
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ECHOSOUNDER OUTPUT 
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